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Abstract 

The reaction ir~p — » nn is investigated within a dynamical coupled-channels model of meson 
production reactions in the nucleon resonance region. The meson baryon channels included are 



c 



ttN, rjN, 7rA, crN, and pN. The non-resonant meson-baryon interactions of the model are derived 
from a set of Lagrangians by using a unitary transformation method. One or two excited nucleon 
states in each of S, P, D, and F partial waves are included to generate the resonant amplitudes. 



> 

£T} ■ Data of tt~p — ► nn reaction from threshold up to a total center-of-mass energy of about 2 GeV are 



satisfactorily reproduced and the roles played by the following nine nucleon resonances are inves- 



^1" 

o 

tigated: <Sn(1535), 5n(1650), Pu(1440), Pu(1710), Pi 3 (1720), Di 3 (1520), Di 3 (1700), I>ib(1675), 

and ^15(1680). The reaction mechanism as well as the predicted rjN scattering length are discussed. 

J3 ■ PACS numbers: 11.80.-m, 11.80Gw, 13.75.-n, 24.10.Eq 



I. INTRODUCTION 

In spite of the quasi extinction of pion beams facilities since about two decades, we are 
witnessing a growing interest in theoretical investigations of pion — nucleon (tcN) inter- 
actions. This is mainly due to the well recognized fact [l|, [2[ that the impressive amount 
of high quality data on electromagnetic meson production reactions from several facilities 
(ELSA, GRAAL, JLab, LEPS, and MAMI) can be used to pin down the underlying reaction 
mechanisms and to study the role and/or properties of intervening baryon resonances only 
when the corresponding hadronic production reactions can also be consistently understood. 
The present work is a prelude to a comprehensive study of the process jp — > rjp, where, 
regardless of the direct production mechanisms considered, a meaningful determination of 
the resonances properties from the rj photoproduction data requires the inclusion of inter- 
mediate and final state meson-nucleon interactions. This latter task is tackled to in the 
present work by analyzing the world data of n~p — > r]n reaction. 

To see the main features of our approach, it is useful to briefly describe here some of the 
recent theoretical works which account for the data of n~p — > rjn reaction. The K-matrix 
coupled-channels approach by Sauermann et al. [3( , included only it N and r/N channels and 
was limited to the Sn partial wave. Such a K-matrix approach was then extended by Green 
and Wycech [4] to include the 7iV channel in a combined analysis of both n~p — > r\n and 
7P — > r\p reactions, and more extensively developed by the Giessen Group [5j, |y, |7| to include 
ttN, rjN, 2nN, ujN, KA, and KT, channels. The approach developed by the Bonn Group 8|, 
[9( is also a K-matrix coupled-channels model supple mented with Regge phenomenology. 



The approach taken by the Zagreb Group 10, llll, ll2l, ll3| has concentrated on performing 



the partial-wave analysis of nN —>■ nN, rjN reactions data. This latter approach is most 
extensively developed by the Virginia Polytechnic Institute - George Washington University 
(SAID) Group [l4|], and is regularly updated. 

In this work, we start with a dynamical coupled-channels model which is based on a 



Hamiltonian formulation [2|] and was applied [15( to analyze nN elastic scattering data. 



This theoretical framework, embodying the nN, rjN, it A, aN, and pN channels, is an 



extension of the approach of Ref. Il6|l and is rather different from the K-matrix models 

nn 

described above, as discussed in detail in Refs. [21. 1151] . Qualitatively speaking, the K-matrix 
approaches, which can be derived [lD from a dynamical formulation by taking the on-shell 



approximation, avoid an explicit treatment of the reaction mechanisms in the short range 
region, where we want to map out the quark-gluon substructure of the excited states (N*) 



of the nucleon. Such a simplification in interpreting the data is a 



so not taken in other 



dynamical approaches such as those developed recently in Refs. |17l . I la . Il9| ] and the earlier 

21] to 



works reviewed in Ref. [l|. Besides the approaches mentioned above, attempts [20, 
introduce subnucleonic degrees of freedom in studying 7i~p — > r\n reaction are also becoming 
available. 

Moreover, combining the dynamical coupled-channels approach and the constituent quark 



model approach [22| to study 



23 



24 



25] -yp —* K + A process proves to be a useful step in 



deepening our understanding of the baryon spectroscopy and search for missing nucleon 
resonances 26] . 



This work follows closely the model (JLMS) developed 15| in a study of tiN elastic 
scattering. The relevant scattering equations are be described in Section II. Section III is 
devoted to model building procedure and evaluation of the data base. In the same Section, 
we present our results for differential and total cross-sections of the process n~p —>■ r]n, in the 
center-of-mass energy range W ^$ 2 GeV, and discuss the main features of the considered 
reaction mechanism. In Section IV the ingredients of the constructed model are used to 
predict the rjN scattering length, as well as the total cross-section for the process r\p — > r\p. 
Summary and conclusions are reported in Section V. 



II. THEORETICAL FRAMEWORK 

A detailed description of the coupled- channels formalism can be found in Refs. [2|, |l5j. We 
outline here the main ingredients which are necessary to understand the procedure followed 
in the present work. 

The meson baryon (MB) transition amplitudes in each partial wave can be written as, 



B 



T M b,m'B'(E) - t MBM , B ,(E) + t MBM , B ,(E) 



(1) 



where, 



MB = ttN, rjN, vrA, pN, aN. 



(2) 



The full amplitudes T mb>M i B i(E) can be directly used to calculate MB — -> M'B' scattering 



observables. The non-resonant amplitude tj^jj M , Br (E) in Eq. ((T]) is denned by the coupled- 
channels equations, 



tMB,M'B'(E) — Vmb,m'B'(E) + 2_, Vmb,m"B"(E) Gm"b"(E) t M „ B ,,^ M , B ,(E) , (3) 

M"B" 

with Gm"b"(E) meson-baryon propagators, and, 

Vmb,m'b'(E) = vmb,m'b' + Z MBM , B ,{E) . (4) 

The interactions vmb,m'b< are derived from tree-level processes by using a unitary transfor- 
mation method. They are energy independent and free of singularities. On the other hand, 
Z MB M i B i (E) is induced by the decays of the unstable particles (A, p, a) and thus contains 



moving singularities due to the ti-kN cuts. As emphasized in Ref. 15[, we neglect that term 
at this stage. 

The second term in the right-hand-side of Eq. ([I]) is the resonant term defined by 

t MB,M>B<{ E )= Yl F MB ^ N *(E)[D(E)]i tj r N *^ M 'B'(E), (5) 

N*,N* 

with the iV* propagator, 

[D-\E)} hJ = {E- M° N: )5 hJ - £y(J3) , (6) 

where M%* is the bare mass of the resonant state N*, and the self-energies are, 

Sjj(-E) = Y ^ n?-*mbGmb{E)T MB ^ N *{E) . (7) 

MB 

The dressed vertex interactions in Eq. (J5J) and Eq. (j7j) are (defining T MB ^ N * = rj v »_ +MB ), 
Tm b^n* (E) = Tmb->n* + 2_^ t MBM , B i{E)GM i B'{E)T M i B i^ N *, (8) 

M'B' 

Fn*^Mb(E) = Tn*-+MB + 2^ ^N*-*M'B'GM'B'{E)t M , B , iMB (E) . (9) 

M'B' 



The parameterization used for T N * iMB is explained in Ref. pj 

The meson-baryon propagators Gmb in the above equations are, 



GmB ^ ^ = F F/Aau- - ( 10 ) 

i5 - E M (k) - £?s(«) + «e 



for the stable particle channels MB = 7riV, r^A^, and, 



GMB(k,E) — — , r , r = y — =r, (11) 



for the unstable particle channels MB = 7rA, pN, aN. The self-energies in Eq. ( JIT]) are 
computed explicitly as defined in Ref. |l5 |. 

To solve the coupled-channels equations, Eq. ([3]), we need to regularize the matrix ele- 
ments of vmb,m'b'- We include at each meson-baryon-baryon vertex a form factor of the 
following form: 



F(k,A) 



k 2 



(12) 



k 2 + A 2 . 

with k being the meson momentum. For the meson-meson-meson vertex of v f , the form 
factor in Eq. (TT2l) is also used with k being the momentum of the exchanged meson. For 
the contact term v c , we regularize it by F(k,A)F(k',A'). Here we follow Ref. [15J and use 
the parameter values determined there for all non-resonant terms except the ones explicitly 
mentioned in the following Sections. 

With the non-resonant amplitudes generated from solving Eq. ([3]), the resonant amplitude 
tf IB M , B , in Eq. ([5]) then depends on the bare mass Mj^* and the bare N* — > MB vertex 
functions. The vertex functions are parameterized in the following way, 

( 2 + L ) r , iL 



N*,MB(LS) 



(k) 



1 



1 



(27t)3/ 2 ^ 



G/v* 



N*,MB(LS) 



A 2 



A 2 



+ {k- k R y 



k 



m T , 



.(13) 



where L and S are the orbital angular momentum and the total spin of the MB system, 
respectively. Cn*,mb(ls) measure the meson- nucleon-iV* coupling strength for a specific LS 
combination of the MB system and are taken as free parameters, and kn are parameters 
fixed by the 7riV — > nN analysis in Ref. [l5j. The above parameterization accounts for the 
threshold k L dependence and the right power (2 + L) such that the integration for calculating 
the dressed vertex Eqs. (IE|) and ([9]) is finite. 



III. RESULTS AND DISCUSSION 



The world data base for the process under investigation embodies 1508 differential and 



98 total cross-sections 



27 



28 



29 



30 



31 



32 



33 



M 



35| for 1.47 < W < 2.85 GeV. 



However, those data presented in 12 papers, thesis, and reports have been obtained mainly 
between 1964 and 1980, except for recent results from the Brookhaven National Laboratory 
and using the Crystal Ball detector by Morrison 27j] and Prakhov et al. 28] . The quality of 
data obtained before 1980 has been discussed by Clajus and Nefkens 36] and emphasized 



by George Washington University [14j, Zagreb [10J and Giessen p, [Z| Groups, underlying 
inconsistencies among different data sets, because of experimental shortcomings and the 
underestimate of systematic uncertainties. This uncomfortable situation has leaded various 

n 

authors to use a reduced data base. For example, the GWU Group [14J] includes in the data 



base 257 data points, mainly from differential cross-section measurements 
but also about 50 data points for total cross-sections 



27 



29 



27. 


28. 


29 



30J, 



30. 3. r 



In the present work, we concentrate only on the differential cross-sections for W < 2 
GeV, as summarized in Table [H The number of data points included in the fitted data base 
in this work (294) will be discussed in Sec. IIII Al 

One of the delicate points in dealing with those data is related to the systematic un- 
certainties (S sys ). For the most recent data by Prakhov et al. 28( those uncertainties are 
given clearly by the authors (6%). For old data, we have mainly followed the general trend 
suggested in Ref . [36J , as summarized in the last column of Table fl] Deinet et al. [29j report 
two sources of systematic uncertainty: 7% and 9%, to be added up quadratically, giving 



Jsys 



11.4%. For Richards et al. [30], we have used 5 sys = 10% for the lowest energies, 



11% to 14% 
Brown et al. 



or other ones, as given in the original paper. For Deben 



32 



ram et al. 31[ and 



, we have followed the conclusion of the Zagreb Group 10, |l2j . In the case 



of Brown et al. 32], we also have lowered the momentum by 4%, in lines with Ref. [36j. 

Total cross-section data has not been included in our fits due to the following reasons: 
i) differential cross-sections are measured by various collaborations in significantly different 
angular ranges with respect to extreme ones (see second column in Table HI), ii) there is 
no commonly agreed upon procedure to extract total cross-sections from measured angular 
distributions, iii) model predictions for extreme angles do not in general agree with each 
other. 



A. Fitting procedure 



As mentioned above, in Ref. 



15| the ttN — > nN reaction was studied within a coupled- 



channels formalism, with multi-step processes embodying ttN, rjN, 7rA, aN, and pN states. 

In that work 175 adjustable parameters were introduced to fit amplitudes produced by 

the SAID Group, fitting more than 10000 data points. About 30 of those parameters are 

particularly relevant to the coupled-channels mechanisms for the ir~p — > rjn reaction. Ac- 



Ref. Angular range P w W N^ p N^ p used in the 5 sys 

(Degrees) (GeV/c) (GeV) the present work 



Prakhov et al. [28] 23 - 157 0.687 - 0.747 1.49 - 1.52 84 

Deinet et al. [29] 32 - 123 0.718 - 1.050 1.51 - 1.70 83 

Richards et al. [30] 26 - 154 0.718 - 1.433 1.51 - 1.90 70 

Debenham et al. [31] 162 - 172 0.697 - 0.999 1.49 - 1.67 111 

Brown et al. [32] 18 - 160 0.724 - 2.724 1.51 - 2.45 379 

Morrison [27] 46 - 134 0.701 - 0.747 1.50 - 1.52 34 

Crouch et al. [33] 14 - 167 1.395 - 3.839 1.88 - 2.85 731 

Feltesse et al. [34] 20 - 160 0.757 1.53 16 



(U 


D7o 


80 


11% 


66 


10% to 14% 


27 


10% + 0.02 mb 


51 


10% or 0.01 mb 



TABLE I: Summary of differential cross-section data for the reaction ir~p — > ijn. Data sets in- 
vestigated in the present work are singled out in the last two columns, where the number of data 
points (Ndp) per data set used in the fitting procedure is given. 

cordingly, in the present work we use that reduced set of adjustable parameters (see Table ILT1) 



and fix the others to their values as determined in Ref. 



15) (cf. Tables III to VII in that 



paper). A total of 294 data points are fitted in the present work (see Tabled]). Here, we 

from the fitted 



wish to make a comment on the exclusion of a few data points in Ref. 
data-base. Actually, as mentioned above, on the one hand, recent data by Prakhov et al. 28] 
bear much smaller errors than other data, and on the other hand, the data-base suffers from 
some inconsistencies. One of the consequences of this situation is that a few data points 
introduce large x 2 s (around 10 or more) thus reducing significantly the efficiency of the 
minimization procedure. The excluded points concern mainly the two lowest energy sets of 



Ref. 



28 



In the following we present our results for two models, as well as those obtained using 
the parameters reported in Ref. [15J (see Table ILT1 and Figs. [EH2]). 

Here, in lines with Ref. [151 ] . the following nucleon resonances (N*) are considered: 
Sn(1535), Su(1650), Pn(1440), Pn(1710), Pi 3 (1720), £>i 3 (1520), Di 3 (1700), Di 6 (1675), 
and F 15 (1680). 

The adjustable parameters for non-resonant terms are the r/NN coupling constant f V NN 
and cut-off A^iv- For resonant terms the parameters are as follows: A^*s bare-masses Mq , 



rjNN* coupling strengths C v nn* , and cut-offs A 



r\ N N* ■ 



Category 



Parameter Model A Model B Ref. [15] 



Non-resonant r)N parameters: 



Bare mass Mq 



* dp 



frjNN 


4.9936 


4.9999 


3.8892 


A-rjNN 


592.11 


591.91 


623.56 


5n(1535) 


1809 


1808 


1800 


5n(1650) 


1901 


1861 


1880 


Pii(1440) 


1775 


1784 


1763 


Pn(1710) 


2019 


2057 


2037 


Pis(1720) 


1726 


1691 


1711 


£>is(1520) 


1918 


1919 


1899 


Dl 3 (1700) 


1971 


1968 


1988 


£>i 5 (1675) 


1878 


1878 


1898 


Fi 5 (1680) 


2207 


2207 


2187 



<^N*^MB(LS) : 












C»7iVSii(1535) 


8.4269 


7.8344 


9.1000 




CpNS-n (1535) 


2.0280 


-0.4935 


2.028 




C»7iVSii(1650) 


2.0487 


-0.4221 


0.6000 




CpiVSn(1650) 


-9.5179 


2.0000 


-9.5179 




C7jAfPn(1440, 


1.6321 


1.6298 


2.6210 




C>7iVPll(1710, 


2.4925 


2.4994 


3.6611 




C 7) JVP 13 (1720, 


2.4474 


2.4997 


-0.9992 




C!7JVD 13 (1520 


) 0.4440 


0.4267 


-0.0174 




C> 7 iVI3 1 3(1700 


, -1.8985 


-0.6463 


0.3570 




C V ND 15 (1675 


) 0.2456 


0.3437 


-0.0959 




C^jvf 15 (1680; 


-0.0446 


-0.0265 


0.0000 


AjV*->MB(LS) : 












^»;]VSii(1535) 


779.38 


799.90 


598.97 




ApJVSn(1535) 


1999.8 


670.89 


1999.8 






1893.8 


955.8 


1893.8 




^»;]VSii(1650) 


500.07 


1999.70 


500.02 




ApJVSn(1650) 


796.83 


2000.00 


796.83 




At?IVPii(1440) 


1766.80 


1757.40 


1654.85 




A,, NP lt (1710, 


500.08 


500.00 


897.84 




A,,jVP 13 (1720) 


631.90 


649.11 


500.23 




^riND 13 (1520 


, 500.20 


500.01 


1918.20 




A,,jV.D 13 (1700 


) 540.55 


763.13 


678.41 




^t?]VDi 5 (1675 


) 507.64 


500.00 


1554.00 




A,, ]VF 15 (1680) 


811.72 


1073.80 


655.87 



2.03 



1.94 



6.94 



TABLE II: Parameters for models A and B determined in this work. The last column gives the 

n 

values determined in Ref. 



>i nit 



Model A is obtained by fitting the data base and those 29 adjustable parameters (see 
column 3 in Table [III)- Model B, for reasons explained below, has five more adjustable 
parameters, namely, the coupling constants and cut-offs of pNSn, with Sn = S , n(1535) 
and 511(1650) for [LS] = [0, 1/2], as well as the cut off A pN s n (i535) for [ LS ] = % 3/21. 



Finally, for comparisons, we reproduce in Table HT1 the relevant values reported in Ref. 15]. 
As mentioned above, in that latter work, adjustable parameters are determined via the 
ttN — > nN channels, and for the n~p — > rjn the data base embodied only a few total cross- 



section data from Refs. 



28 



32( | . Notice that the five pNSn parameters in model A (shown 



in italics in Table HT1) were not treated as adjustable parameters, and hence, are identical to 
those of Ref. [13]. 



B. Differential and total cross-sections for the process it p — > r\n 

In Figures Q] and [2] we compare the results of the models A and B with the differential 
cross-section data, for which the reduced x 2 s per data point are 2.03 and 1.94, respectively. 
Those numbers compare well with the GWU [l_4| reduced \ 2 = 2.44. In the same figures, 



we also show results obtained by using the parameters of Ref. 15j, which gives \ 2 = 6.94. 

Before discussing different curves in comparison with data, we wish to emphasize the 
difference between models A and B. Once the model A obtained, we investigated the 
importance of various parameters and found that by switching off the p coupling to the 
Sii(1535), the x 2 increases by roughly a factor of 3. Within the investigated reaction, such 
a high sensitivity to the pNSn seems unrealistic. To cure that behavior, we refitted the data 
by allowing those coupling constants to vary in the range of ±0.5 for S , n(1535) and ±2 for 
Sii(1650), instead of ±10. The model B is then obtained, where that effect is significantly 
reduced. Comparing the two models in Figs. [T] and [2J we observe that they differ from each 
other by less than the statistic uncertainties of the data, corroborating that the 7r~p — > r]n 
reaction is not a proper channel to pin down those couplings. Better constraints on those 



couplings can be obtained by investigating the %N — > -miN process [37 1. 

Models A and B show reasonable agreements with Prakhov et al. [28] data, except at the 
lowest energy (Fig. [ p. We will come back to that point later. At four common energies, 
27j are also depicted. That latter data set, not included in our fitting 



data from Morrison 



procedure, shows systematically smaller cross-sections compared to Ref. 



data. 




-1 -0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 -1 -0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 

cos cos 



FIG. 1: Differential cross-section for the reaction it p — > rjn. The curves correspond to models 
A (dashed) and B (full) from the present work. The dash-dotted curves are obtained by using 



. m 



la. Data are from Prakhov et al. 281 (empty circles), Deinet et al. 29] 



the parameters in Ref. 

(cross), Richards et al. [30(] (empty squares), Debenham et al. [31] (up-triangles), and Morrison 

(diamonds). Data uncertainties depicted are only statistic ones. 



27] 



Prakhov et al. 



281 ] data set at W = 1.507 GeV is of special interest, since there are also 



data from three other measurements. Results from Deinet et al. 29| are compatible with 
Prakhov et al. data, though with larger uncertainties (which become even more sizeable 
at W — 1.525 GeV). Richards et al. [30] data show deviations from Prakhov et al. ones, 



especially at most backward angles. Finally, copious data from Debenham et al. 



are 



unfortunately limited to extreme backward angles and appear to be rather consistent with 



10 




-1 -0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 -1 -0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 
COS COS 



FIG. 2: Differential cross-section for the reaction n p — ► r\n. Curves anc_ 
Additional data are from Brown et al. [32] (right-triangles) and Crouch et al. 



data as in Fig. CD . 
331 ] (down-triangles). 



other data only up to V7 ~ 1.55 GeV. This trend is confirmed in Fig. [2], where models A 
and B reproduce correctly results from Deinet et al. 29] and Richards et al. [30| . Both sets 
of data come out fairly compatible with measurement from Brown et al. [32] at W = 1.699, 
1.729, and 1.805 GeV. Models A and B show acceptable agreements with those data, except 
at backward angles, where the model / experiment discrepancies get reduced when energy 
increases and suitable agreement is observed at W = 1.897 GeV. At the three remaining 
depicted energies (W = 1.871, 1.948, and 2.003 GeV) our models reproduce the general 
trend of Brown et al. data. At those energies, data from Crouch et al. 33j, not included 
in our data base, are also shown. The two data sets are not consistent. Given the known 



problems 36) with Brown's data, we made also attempts to fit the data base, within model 



11 



A, by replacing the Brown et al. data by those of Crouch et al. at W = 1.879 and 1.915 
GeV. However, we observed a significant increase of x 2 which goes from 2.03 to 4.12, and 
with very undesirable effects in the Crystal Ball energy range. 
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FIG. 3: Total cross-section for the reaction ir~p — > rjn. Curves are from Ref. [15( (dash-dotted), 
model A (dashed), model B (full), and the background contributions (dotted) in model B. Data 
are as in Figs. [Qand[2j 



In Figures [T] and [2j results using the parameters in Ref. [ToJ are also shown. At lowest 
energies, that model overestimates the data. At higher energies, it shows significant devia- 
tions, first at backward angles and then at forward angles. Above W « 1.8 GeV it tends 
to miss the data. 

Finally, as mentioned above, we did not include the extracted total cross-section data 
in our data-base. In Fig. [3l we show the postdictions of our models A and B, as well as 
results of the Ref. 15j , and compare them with the data. Both models A and B reproduce 
correctly the data, except for those by Crouch et al. 33j, for which the differential cross- 
sections turn out to be significantly smaller than other data, as shown in Fig. [21 Moreover, 
the background contributions show a smooth behavior and are small with respect to the full 
model results, except around the minimum of the total cross-section, where resonant terms 
produce highly destructive interferences. 



To end this Section, we wish to emphasize that the results of Ref. 15] are extended to 
the process ir^p — > rjn and two models are obtained, reproducing equally well the general 
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FIG. 4: Comparisons between the results from Ref. [15|] (dashed curves) and from the model B 
(dotted curves) for ttN — ► X, 7riV processes. Left panel: Predicted total cross-section for the 
reactions n + p — > X (upper set) and n + p — » n + p (lower set). Right panel: predicted total cross- 
section for the reactions n~p — ► X (upper set) ir~p — > ir~p + 7r°n (lower set); Data in both panels 



are from Refs. 



38 



39] 



features of a heterogeneous data-base. This new set of parameters, particularly relevant 
to the investigated process, does not spoil the excellent results obtained in Ref. [15[, and 
devoted mainly to the irN — > nN observables. In order to illustrate this latter point, results 
from Ref. 15[ and our model B are shown in Fig. HI where Figs. [13] and [14] of Ref. [15| 
have been complemented with the predictions of the model B. For the processes with n + p 
initial state (left panel in Fig. H|), results from the two models overlap with each other. For 
reactions involving ir~p initial states (right panel in Fig. H]), model B gives deeper minima 
around W « 1.4 GeV than those reported in Ref. 15| . with the largest discrepancy between 
the two curves being less than 20%. Comparing the partial- wave amplitudes of those models, 
the main differences appear in the P\\- and Pi3-waves. 



C. Main features of the tt p — > r\n reaction mechanism 

In order to shed insights to the main ingredients of the reaction mechanism, we concen- 
trate on model B. Starting from that model, and without further minimizations, we have 
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checked the variations of the x 2 by switching off the nine resonances one by one. The results 
are reported in Table IIHI 

5n(1535) £n (1650) Pu(1440) Pn(1710) Pi 3 (1720) £>i 3 (1520) £>i 3 (1700) Di 6 (1675) Fi 5 (1680) 
X 2 48.86 2.62 3.55 2.37 2.77 2.23 1.93 2.10 2.47 

TABLE III: Reduced x 2 P er data point for model B with one resonance switched off (the reduced 
X 2 for the full model B being 1.94). 
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FIG. 5: Differential cross-section for the reaction 7r~p — ► rjn. The full curves correspond to model B 
and the dotted ones to the non-resonant terms contributions. The other curves have been obtained 
by removing one resonance from that model; the removed resonances are: S r n(1535) (dashed), 
Pn(1440) (dash-dotted), and Pi 3 (1720) (dash-dot-dotted). Data are as in Fig. [TJ 



14 




-1 -0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 -1 -0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 
COS COS 



FIG. 6: Differential cross-section for the reaction n p — > rjn. The curves are is in Fig. Data are 
as in Fig. [2j 

As expected, the process is dominated by the Sii(1535) resonance. There are however 
two other resonances playing non-negligible roles, namely, Pn(1440) and Pi3(1720). Figures 
E] and M show that the importance of those resonances depends on both angle and energy, 
^he Sii(1535) resonance produces more than 80% of the cross-section for the Prakhov et 



al. 



28j data. Its importance decreases with energy, especially at backward angles, without 



vanishing out. The effect of the Pn(1440) resonance becomes visible roughly in the energy 
range 1.525 < W < 1.8 GeV, with a destructive behavior at most forward angles. Finally, 
the Pi3(1720) appears, in the forward hemisphere, around W ~ 1.6 GeV, with the highest 
contributions at W ~ 1.73 GeV and its effect remains constructive. 

Although the effect of the _Di3(1520) on the y 2 is small, it is required to produce the right 
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curvature of the curves at low energies. 

In conclusion, model B turns out to describe in a satisfactory manner the data set and 
embodies a simple reaction mechanism. In the following Section we use hence that model 
for further investigations of the r]N system. 

IV. PREDICTIONS FOR THE r]N SCATTERING LENGTH AND THE r/p -> rjp 
TOTAL CROSS-SECTION 

The rjN scattering amplitude in terms of the t-matrix is given by the following relation: 



/(*) 



-TV 



yk 2 + m 2 N Jk 2 + 7?7,2 



k 2 + m 2 N + Jk 2 + ml 



t nN (k,k). 



(14) 



Then, the scattering length reads, 



a vN = lim k -Mif(k). 



(15) 



Fig. [7] shows the real and imaginary parts of the function f(k), for model B, and leads 
to the following value for the scattering length: 



a vN = (0.30 + z0.18) fm. 



(16) 



Efforts since about two decades to determine the i]N scattering length have recently been 



reviewed by several authors 40 



41 



42J . A lower limit on t 



re imaginary part, derived from 



27J leads [42j to 0.172±0.009 fm. 



the optical theorem and taking into account the recent data 

Combining results quoted in those review papers, the present knowledge of the imaginary 

part is: 



0.17 < Im a vN < 0.49 fm, 



and our value comes out to be within that range. 



(17) 



For the real part of the scattering length the estimates in the literature give 41] 



0.27 < Tie a vN < 1.0 fm. 



(18) 



The value extracted in the present work, within model B, is close to the lower limit. Our 
value is compatible with those obtained via chiral effective Lagrangians 43|, most recent 
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FIG. 7: r/N scattering amplitude f(k) as a function of cm. momentum, within model B. 



solution (G380) from energy-dependent partial-wave analysis [42JL of e 



and n p — > rjn scattering data, as wel 



on chiral perturbation approaches 46 



as with older findings [44j, |45[. Investigations based 



astic tv p, it p — > n°n, 



47| lead to smaller values around 0.2 fm. Finally, 
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48] or K-matrix 49| produce larger 



coupled-channels calculation within T-matrix 10, 
values, 0.5 < IZe a vN < 1.0 fm. 

Besides the process ttN — > r]N, the 77 production using proton or deuteron beams have 
also been investigated using various sets of r]N scattering lengths reported in the literature, 
as summarized below: 

a) pn — > rjd near threshold data [50j has been studied within a two-step model 45j, em- 



bodying meson-exchange and final-state r]N interactions, and favors small scattering len gth : 



a>r)N = 0.29 + i0.36 fm. A microscopic three-body approach in its non-relativistic version 51] 
reached the conclusion that the data are well reproduced using the results from Ref. [46], 
o-tjN = 0.42 + i0.34 fm. The relativistic version of that approach [52] shows the importance 
of initial- and final-state treatments, emphasized also by Julich Group 53), leading to a 
reduced selectivity on the sets used for the scattering length 
b) pp 



54) and pn — > pnrj [50] data, as well as the above mentioned data have 



pprj 



recently been studied within an effective Lagrangian model 55j resulting in a reasonable 
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account of data for a 



i)N 



0.51 + i0.26 fm. 



c) r\ production in proton-deuteron collisions are being studied |56|, but at the present 
time those investigations do not allow refinements in determining the i]N scattering length. 
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FIG. 8: Total cross-section for the reaction rjp — ► r/p as a function of total cm. energy, within 
model B. 

Findings of various approaches with respect to the r]N scattering length, summarized 
above, lead then to the ranges for real and imaginary parts as reported in Eqs. (1171) and 
( TT8"j) . Our value for the real part being close to the lower bound, excludes the existence of 
bound //-nucleus states. 

To end this Section, we show our prediction for the rjN elastic scattering total cross- 
section (Fig. [8]). The background contributions (dashed curve) turns out to be small and 
smoothly varying. This latter contribution completed by that of the Sii(1535) (dash-dot- 
dotted curve) account for a significant portion of the total cross-section predicted by the full 
model B (full curve). 

V. SUMMARY AND CONCLUSIONS 

A dynamical coupled-channels formalism is used to study the the process n~p — > rjn 
in the total center-of-mass energy range W ^ 2 GeV. The formalism embodies, besides 
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non- resonant terms, five intermediate meson- nucleon states, namely, nN, rjN, it A, aN, and 
pN. 

Within this phenomenological approach, 34 adjustable parameters are introduced, two 
of them for the non-resonant mechanisms and the others for the nine nucleon resonances 
retained in the model search, namely, 5 U (1535), Sn(1650), P u (1440), Pn(1710), Pi 3 (1720), 
£>i3(1520), Di 3 (1700), -Di 5 (1675), and Pis(1680). That set of resonances corresponds to all 
known 3 and 4 star resonances relevant to the energy range investigated here. 

In order to determine the parameters and build a model, a data set including 294 measured 
differential cross-sections, coming from five collaborations, are fitted. The selection of data 
points allows to suppress the manifestations of inconsistencies among available data sets. 
Our model B reproduces satisfactorily the data, with a reduced x 2 — 1-94. A detailed 
study of the reaction mechanism within model B allows to establish a hierarchy in the 
roles played by nucleon resonances. Actually, the dominant resonance turns out to be the 
Sii(1535). The other resonances affecting the x 2 by more than 20% when switched off, 
are by decreasing importance: P u (1440), Pi 3 (1720), 5 u (1650),Pi B (1680), Pn(1710), and 
.Di3(1520). Contributions from Pi 3 (1700) and Di 5 (1675) are found to be negligible. 

Model B is used to extract the i]N scattering length, which comes out to be: a V N = 
(0.30 +i 0.18) fm. Both the real and imaginary parts of that quantity are within the ranges 
determined from other works. 

To improve our knowledge on the n~p — > rjn, and consequently on the r]N system, 
further measurements including polarized target asymmetry, are highly desirable. Such 
experimental results are awaited for thanks to present and/or forthcoming pion beams in 



the following labs: GSI [57], ITEP [58J, Fermi Lab |59j, and JPARC [60|. 



61]. 



Finally, to take advantage of copious rj electromagnetic production data, the obtained 
model B appears reliable enough for our in progress investigation of_the jp — > r]p reaction 
within a coupled- channels approach and a constituent quark model 
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